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We discuss the possible existence of a "spin-gap" phase in the low-doping regime of strongly- 
correlated two-dimensional electrons within the gauge field description of the t — J model. The 
spin-gap phase was recently shown by Ubbens and Lee to be destroyed by gauge field quantum 
fluctuations for a single-layer 2D system in the absence of disorder and for a full gap. We show that 
the same conclusion applies both in the dirty limit and for the case of a gapless spinon condensate. 

PACS Nos. 74.30E; 64.70D 



I. INTRODUCTION 



Strong electronic correlations have been described as essential to the understanding of the quasi-two-dimensional 
high-T c superconductors. 1 Spin-charge separation is a possible manifestation of strong interactions which has been 
proposed as basic to the description of the low-energy excitations in the 2D interacting electron system. 2 For strong 
interaction, the t — J model 3 version of the Hubbard model is often used as a starting point. The system is studied 
when it is doped with x holes per site away from half- filling (thus 1 — x electrons per site). Configurations in this 
model are constrained to exclude those where the number of electrons on a given site can exceed unity. The slave 
\Q boson method 4 is a way to deal with this strong-correlation constraint and it leads naturally to spin-charge separation 
ON via the introduction of x charged bosonic holons and 1 — x chargeless fermionic spinons. A mean-field approach yields 
the phase diagram sketched in Fig. 1. It contains a number of phases including a "strange metal phase" and phases 
in which either the spinons or the holons (or both) are condensed. 5 It has been suggested 6 that the properties of 
the normal state of high-T c superconductors can be well-described by the strange metal phase which is essentially a 
uniform RVB state. 1 . 




FIG. 1. Mean field phase diagram of the t — J model in the temperature (T)-doping (x) plane, containing the "strange" 
metal phase (I), Fermi liquid phase (II), superconductor phase (III), and spin gap phase (IV). 

A convenient way to discuss the phases is to consider the Ioffe-Larkin expression 7 for the response function II of 
the total system to an external electromagnetic field: 

where ILf and II;, are the current-current response functions of the fermion (spinon) and boson (holon) subsystems. 
This equation is a consequence of the strong-correlation constraints. The response of the total system is dominated 
by the smaller of ILf and TIj, . Accordingly, the schematic phase diagram shown in Fig. 1 has four parts. When both 
fluids are normal, i.e. when the spinons are a normal Fermi gas and the holons are not Bose-condensed, as initially 
studied in Ref. 4, we have a strange metal, with a Fermi surface for 1 — x electrons which obeys Luttinger's theorem, 
but with transport properties determined by the doping x. Thus there is a positive Hall constant and a resistivity 
proportional to \ jx. The temperature dependence of the resistivity, p oc T 13 with (3 near unity, is determined by the 
neglect of the holon diamagnetic susceptibility in comparison with the spinon diamagnetic susceptibility which enters 
the denominator of the n response functions. 
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The mean-field theory of the system accounts for the local constraints only globally. The corrections involve the 
fluctuations of a gauge-field which mediates the conservation laws of the bosons and fermions which are implied by the 
constraint. The reason for the unusual resistivity law is that the less-conducting particles, the bosons, are inelastically 
scattered by low-frequency long-wavelength gauge fluctuations. 

The retarded gauge field propagator D(q,u>) is given by: 



d- 1 = n f + n b 

n/,6 = Xf,bq 2 - iw<jf,b(q, w). 

The spinon (fermion) and holon (boson) susceptibilities are 

1 

Xf 
Xb 



(1.2) 
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In the pure case, the spinon and holon conductivities are 

2n f 



<r f (q,u)) 
Vb(q,v) 



m f qv f 

m b qv b ' 



oj < qvf 
w < qv b = q-\/2T/m b , 



(1.4) 



where v b = yjll ' jm b is the thermal velocity of bosons. In the above equations, the fermion and boson areal number 
densities and masses appear. They are given in terms of the energy parameters of the t — J model and the doping x 
away from half-filling as 



{n t lm t ) = {k 2 F /2-Km t )~2{\-x)J 

(n b /m b ) ~ 2xt, 



(1.5) 



where kp is the Fermi momementum of spinons. The conductivities are zero for w > qvf tb . In Eqs. (1.2), we have 
assumed that the temperature is sufficiently high that the bosons are far from condensation and we have treated 
their distribution function classically. For non-interacting bosons, this would require that 2xt/T < 1. However, the 
interactions extend this behavior to a wider temperature range. 8 

In the strange metal phase, the fermionic contribution dominates the gauge fluctuations and we have 



2(n f / mf )[-i(w/qv f ) + (q 2 /12k 2 )] 



(1.6) 



Ubbens and Lee 9 have recently argued that in the pure system the spin-gap phase (region IV in Fig. 1) is an artifact 
of the mean field approximation. According to them, the spin-gap is destroyed by the quantum fluctuations of the 
gauge field. Their argument goes as follows: The possible condensation of the spinons into the spin-gap phase is 
like that of the superconducting transition. As in BCS theory, one can compute the negative contribution F conc i(A) 
of the condensation to the total free energy as a function of the spinon pairing amplitude A. However, the free 
energy has a contribution from the gauge field as well 10 because when a spinon gap opens, a mass is induced in the 
fermionic contribution to the gauge fluctuation spectrum and the low-lying gauge field fluctuation contribution to the 
free energy is lost. This gives a positive contribution to the free energy. It is 9 ' 10 



Fgauge(A) F gauge (Q) 



J (27T) 2 J 



du> coth 



( — ) arctan 
\2TJ 



lmD~ 1 (q, u, A) 



ReD~ 1 (q, w, A) 



-fgatige(O). 



(1.7) 



It is necessary to examine the effect of A on _D _1 = T[f + T[ b . In lowest order, only is affected. As in 
superconductivity, the massless of the normal phase [see Eq. (1.2)] acquires a gap proportional to the superfluid 
density (Meissner effect). We have n^(A) = n s f/rrif + 0{w,q 2 ). Here n s f = nf at T = in the clean (London) 
limit and n s f = UfArf in the dirty limit (At^ < 1), where Tf is the spinon elastic transport lifetime (see Ref. 11, for 
example). In fact, ImILf(q,u!) and Tlel[f(q,ui) become complicated functions of w, especially when the gap function 
is anisotropic. The main effect on the low frequency part of is that Imn^ = for w < 2 A at T = 0. With neglect 
of possible spinon gap anisotropy, an estimate of F gauge (A) — F gauge (0) in the clean limit is given by the contribution 
to F gauge (Q) of the gauge-field modes with w < 2 A which are missing in the condensed phase 9 : 
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/an /\ f , w w n /" d 2 q ( ntw \ 

F gauge (A) - F gauge (0) ~ ^ dw coth( — ) y ^-^ arctan (^^^ J 



x 2 / 3 

\ 2/3 

) A 5 / 3 , (1.8) 

where x is the coefficient of q 2 in Re_D _1 , x = (l/127rray) + (n5/12rajT). Since F gauge (A) ~ A 5 / 3 is non-analytic in 
A 2 , the spinon pairing transition must be actually be first order. More importantly, F gauge (A) — F gauge (0) is always 
much larger than |_F co)l( ;(A)| for any A, as long as x remains small. This implies that the opening of a gap does not 
happen for T > Tbe = 2'irni > /mi > . For completeness, we have reproduced here the argument of Ubbens and Lee 9 but 
we included the holon contribution to Re_D -1 in Eq. (1.8). This reduces very slightly the impact of F gauge but does 
not change the conclusion. 

One might wonder if the vanishing of the mean-field spin gap is a property of the pure system alone, since low- 
frequency gauge fluctuations are significantly modified in the dirty limit. The conductivities will be finite as q — > 0, 
contrary to the pure case. As a result, the zero-point fluctuations should be significantly reduced in the long wavelength 
limit. The most favorable situation for the preservation of a mean-field spin pseudogap would be a gapless situation, 
whereby a finite density of excited states continues to exist down to zero energy. 

The purpose of this paper is to examine this question. Our conclusion is that the mean-field spin gap phase is 
also suppressed in the dirty limit of the spinon superfluid, as well as in the gapless case. Although there remain low 
frequency gauge fluctuations associated with the holon liquid, for T > Tbe, the suppression of the low frequency 
modes in the normal phase proves too costly, in comparison with the gain in spinon gap condensation energy, to allow 
the condensation to happen. 



II. DIRTY LIMIT 
A. Normal phase 



In the normal phase, disorder produces a non-zero long- wavelength conductivity for each component of the plasma. 
In a Fermi liquid, electron-impurity scattering gives a conductivity nyr e /my with an elastic scattering rate T e -1 = 
2irniNo\V\ 2 . Here, V is the Fourier coefficient of the weak short range electron-impurity potential, n, is the density 
of impurities and No is the one-spin density of states. The effect of impurity scattering in the case of the strongly 
correlated electron gas has been treated extensively by Ioffe and Kotliar 12 (IK) in the context of the t — J model. There 
also exist discussions for the Hubbard model. 13 Following IK, we express the electron scattering rate in terms of the 
variance w of the random potential: T e -1 = 2irNow 2 . It was pointed out by IK that in the strongly-correlated electron 
gas, impurity scattering results in quantitatively different transport lifetimes for the fermions and the bosons. The 
interplay between the non-uniformity in space of the random potential and the strong-correlation constraints creates 
non-uniform chemical potentials (Lagrange multipliers) for the spinons and holons. The holons see an unscreened 
random potential, while the effective spinon-impurity potential is screened by the backflow of the holons; it becomes 
smooth and only forward scattering is appreciable. Consequently, the spinon transport lifetime is much longer than 
the boson one: 

247rf 2 

Tf ~ Jw 2 x*l 2 ' ( ' 

while t f, ~ 2t/w 2 , so that Tf /ti ~ 12trx~ 3 / 2 (t/ J) » 1. As a result, as emphasized by IK, the resistivity is determined 
by the holons. 

The fermion conductivity <jf = nfTf/nif is determined by the above value of Tf as long as qvfTf < 1. Otherwise, 
the pure value given in Eq. (1.4) obtains. In fact, the fermion conductivity departs from its value in the pure case 
only in a small region of reciprocal space, i.e. for q/kf < 10 -2 to 10 -3 , if we choose w ~ t and x ~ .1 — .2. On the 
contrary, the boson conductivity is governed by impurity scattering in the whole reciprocal space. 

We conclude that the gauge field propagator D(q,u>) = (q, u>) + ILb(q, w) is practically determined entirely by the 
spinon polarisability, while the physical conductivity is dominated by the holon fluid: 

- - GfGh ~ a b = 4x(t/w) 2 . (2.2) 



<7f + &b 



3 



Then the free energy cost due to a mean field spin gap A has the form originally given in Ref. 9 and quoted here in 
Eq. (1.8): 



-F n = 38.37\T(0)A 2 
where N(0) is the one-spin density of states, rrif /2w. 



(1 - x)J 



A 



1/3 



(2-3) 



B. The spin gap phase 



In the spin gap phase, in the dirty case, a free energy contribution from gauge fluctuations is recovered, the fermion 
polarisability is that of a dirty superfluid, i.e. Imllj: = 0, Rellj: = (nyA-ry/ray) + 0(w), valid if Arf < 1. Since 
the boson polarisability is unchanged compared to the normal phase as long as T > Tbe, the inverse gauge field 
propagator is now: 



D(q,w) 1 ~ (nf/m f )ATf - iwa b + xq 2 
As in Eq. (1.7), we need the ratio of the imaginary and real parts of _D _1 . We have 



(2.4) 



WXtTf, 



ReD -1 ~ A(l - x)jT f 

xq 



q<q* 

q>q* 



2Ar f (l - x)J] 1/2 



(2.5) 



where we have used Eqs. (1.5). 

Therefore the recovered gauge free energy in the spin gap phase is: 



, w2, r 1 d q f , , f"\ \lmD~ 1 (q, w, A) 

F s + N(0)A 2 2 = / — ^ / dw coth — arctan „ „ \{ , 



(2.6) 



where — 7V(0)A 2 /2 is the contribution to the free energy of the spinon condensation. The susceptibility is still 
dominated by the fermion contribution, x — (1 — s)J/(3fc|). For an example, we take x = .2 and find, in the limit 
A » T: 



F s = -3.067\T(0)A 2 - .57\T(0)A 2 . 
Finally, we total the free energy cost (for x = .2) to condense the spinons: 

F s -F n ~ -35.57V(0)A 2 [.100- (J/A) 1/3 ] 



(2.7) 



(2.8) 



Thus, in the dirty limit, the mean field spin gap can only exist if A >> J. It is straightforward to check that the 
transition can only be first order. 



III. GAPLESS CASE 



Gapless superconductivity can arise in a variety of ways. In the context of s- wave superconductivity, scattering from 
magnetic impurities causes breaking of singlet pairs and this leads to gapless superconductivity. 11 ' 14 When Ar s < 1 
(t s is the spin-flip scattering rate), the density of states is non-zero at the Fermi level in the superconducting state 
and the dc conductivity has a finite value at q = 0. The superfluid density is reduced and the electron polarisability 
contains a conductivity-damping term iu>a(q, u>) down to zero frequency. When Ar s < 1, <j(q,u>) is non-zero at q = 0, 
and grows linearly with frequency. 

We shall assume that some spin-flip mechanism is at work in the strongly correlated medium, and that the spin-flip 
rate for the spinons is such that the condition Ar s < 1 is satisfied. Here A is the fermion spin gap. Then the spin-gap 
phase behaves as a gapless superfluid. For simplicity, we take Ar s ~ 1. The conductivity in this case 14 can be 
approximated by (uhj„/6A) for w < 6 A and by a„ otherwise. With this choice of a, we obtain the superfluid density 
by applying the /-sum rule: 7rn/2ra = J <j(w)dw, so that 
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s /2m = J (<r ln - cr ls )dw (3.1) 

■/ 

Jo 



. , waw. 
6A 



Therefore, in the spin gap phase, we have 



(this is strictly valid for Ar s < 7r/6) and 



n s QnfAr 
rrif wrrif 



6A 36A 2 



(3.2) 



w < 6A (3.3) 



As usual, in order to assess the stability of the spin gap phase, we compare its energy with that of the normal 
phase. We look at the cost in the gauge field contribution to the free energy which is paid upon entry into the spinon 
condensate phase. It is again given by F Stgauge — F„ tgauge , where s, n stand for the spinon condensate and normal 
phases. Usually, F Stgauge = F„ tgauge for frequencies larger than 2A. In the gapless situation, we have seen that the 
spinon conductivity in the condensed phase is only equal to that in the normal phase for frequencies above 6A. Then 
the gauge field free energy differs between the normal and condensed phases up to a frequency 6A, after which the 
conductivities (hence Imll) in the two phases are equal. Therefore, we have to evaluate the free energy contributions 
up to frequency 6A. In the normal phase, Imll ~ w<jf = w(ny /my)r 4 and Rell ~ Xf<l 2 = (<? 2 /127rray). The spin-flip 
scattering rate is not reduced by screening effects. In the normal state, we have a contribution which differs from that 
of the previous sections since the absence of screening makes the fermion conductivity (now given by nfT s /rrif) exceed 
the Landau damping [given by (rif /rrifVfq)] over a larger portion of the Brillouin zone. For simplicity we take the 
mean free path of order the lattice spacing, implying that the spin-scattering conductivity determines Imll over the 
whole zone. This is an extreme case which favors the spin gap phase. Consequently, we write, in the normal phase, 

Imll = u>T s n f lm f 

Ren q 2 /12irm f ' { ' ' 

To determine the free energy loss from the normal phase, we use this expression in the RHS of Eq. (1.7) and integrate 
the frequency to 6A. The result is (for x = .2 and Ar s = 1) 

F n = -2717V(0)A 2 ^-. (3.5) 

We again recover a gauge field contribution in the condensed phase since the conductivity does not vanish there. The 
calculation of the gauge field free energy in the spinon condensed phase is similar to that of Sec. IIB. In the present 
case however, the fermion conductivity [Eq. (3.3)] is not zero for w < 6A. In fact, for x = .2, it is larger than the 
boson conductivity which was used in Eq. (2.3) as soon as w > .004A, so we use 

Imll = -uht. ~ -I^S. (3. 6 ) 

for all u> up to 6A. 

At the same time the superfluid stiffness adds a constant part to Rell:, 

Ren- -(^At s ) + X9 2 - (3.7) 
7r rrif 

This first term on the RHS exceeds the second over the whole zone. Consequently, in the spinon condensate phase, 
we have 

Imn / w \ 2 .„ „. 

Ren = ^ (6A) • ( 3 - 8 ) 

We use this in Eq. (2.5) and replace the upper limit of the w integration by 6A. The result is 

F., gaug e = -95.67V(0)A 2 (J/A) (3.9) 
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The condensation energy is much reduced in the gapless phase compared to the pure case; with our choice Ar s f ~ 1, 
it is reduced to: 14 

F cond = -MN(0)A 2 (3.10) 

The net free energy cost is 

F s , gaug e + F cond - F n , gauge = - 175tf (0) A 2 [.0003 - (J/ A)]. (3.11) 

Thus, we see once again that the cost in zero point energy exceeds by far the gain in condensation energy since the 
latter is so small in the gapless case. 

IV. CONCLUSION 

The mean field spin gap is destroyed by the gauge fluctuations, in all the cases considered so far: pure case, 9 
dirty limit, gapless phase (this work), in the absence of singular interactions. We have not considered the case of 
an anisotropic spinon condensate, say "d-wave". In the case of a condensate with nodes, there are again gapless 
excitations and a reduction in condensation energy. However, for the case of point nodes, as in the two-dimensional 
d-wave case, the condensation energy remains of order A 2 . Furthermore, as Ubbens and Lee 9 have shown, the presence 
of a non-zero spinon density of states and the fact that for some momentum transfers the gauge field-spinon interaction 
is pair-enhancing rather than pair-breaking does not change their conclusion. Thus, we do not expect a result which 
differs from that of Section III. It might be argued that since impurity scattering is pairbreaking for a gap function 
with nodes, it would effectively destroy a d-wave spinon condensate anyway. However this effect arises when the 
scattering connects two regions of the Fermi surface which have gap functions of opposite sign. However, as discussed 
in Section IIA, the spinon-impurity scattering is predominately of small momentum transfer and therefore most of 
the scattering connects gap regions of the same sign. 

We have not discussed the contribution of longitudinal gauge fluctuations. This might seem dangerous in view of 
the findings of Ref. 10 according to which longitudinal fluctuations have a dramatic effect at high temperature on the 
total entropy of the spinon-holon system; they reduce this entropy by roughly a factor two at low doping. However, 
at low temperature, longitudinal fluctuations are frozen, as in an ordinary charged electron liquid, so that they have 
exponentially small contributions in both the spin gap phase and the strange metal phase. Therefore we do not think 
they change significantly the picture obtained when considering transverse fluctuations alone. 

It is not clear whether the disappearance of the spin gap phase in the gauge theory of the single-layer material 
improves the agreement between that theory and experiment. Ubbens and Lee 9 adopt the point of view of Millis 
and Monien, 15 who argued, on the basis of NMR and neutron diffraction data, that the spin gap phase exists in 
the two-layer compound YBCO (yttrium barium copper oxide) but is absent in the single-layer LSCO (lanthanum 
cuprate). Then they proceed to show that an interlayer pairing mechanism in YBCO restores the spin gap in that 
compound. However, the evidence against the spin gap in LSCO is not completely unambiguous. It has been argued 
that similar suppressions in magnetic, transport and thermodynamic properties is rather generic in all underdoped 
materials 16 ' 17 . It is tempting to ascribe these observations to a spingap mechanism. However, according to the results 
of Ubbens and Lee, as well as of this paper, the slave boson gauge theory cannot account for this suppression of 
inelastic scattering processes in terms of suppressed spin degrees of freedom in the single layer materials. 

While there is certainly disagreement at present, if it turns out that both LSCO and YBCO have quite similar "spin- 
gap" phenomenology, it will be difficult to account for the behavior by different mechanisms for the two compounds. 
Then the inability of the present slave-boson gauge-field approach to describe a spin gap in the normal phase of LSCO 
will be a serious shortcoming of that approach and raise questions about its use to analyze normal state properties of 
YBCO. 5 ' 18 ' 19 

In a recent work, Wen and Lee 20 have introduced a new mean-field theory for the single-layer t — J model. There 
is a staggered flux phase which has a spin gap. Since the gauge field in the approach remains massless in the spin gap 
phase, the latter may turn out to be stable even after consideration of the fluctuations. 
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